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Abstract 


A simulation mode! of Single-mode Single-polarization Vertical Cavity Surface Emitting 
Laser (VCSEL) is developed by which the turn-on event of the VCSEL can be simulated 
by solving the rate equations numerically after introducing the parameters of the VCSEL 
to be simulated in the laser rate equations. The numerical solution of the coupled rate 
equations is done by 4' h order Runge-Kutta method using the software LabVIEW, The 
spontaneous emission noise is simulated by using Langevin noise sources in the 
equations. From the statistical behaviour of the tum-on delay of the zero-biased laser, the 
probability density function(PDF) of the tum-on delay is calculated and plotted for 
different on-states. The turn-on jitter is calculated which is the standard deviation of the 
PDF. Also the relaxation resonance frequency is c alculated for different on-states. Then 
the variation of jitter with relaxation resonance frequency (i.e. the variation of jitter with 
different on-states) is compared with the analytical expression derived from the theory 
and also with the measurement results. Also the effects of different off-states on jitter is 
investigated. Further the effects of different bit patterns and also the effects of bit rates on 
jitter is investigated. 

For the simulation of the VCSEL, the data rate taken for simulation of jitter is 160 
Mbits/s. Also the simulation is extended to 1 Gbits/s. First, the 1010 pattern is applied 
and then the pseudo-random pattern is applied to observe the bit pattern effect. 

The laser characteristic graph obtained from simulation is matching well with the 
experimental graph although the effect of suppressed polarization is neglected in 
Simulation. Also, the variation of jitter with relaxation resonance frequency (f r ) is in 
good agreement with the measurement results as well as with the theory. When the 
relaxation resonance frequency exceeds 3 GHz, the variation of jitter with f r is not 
matching with the experimental results as the laser is multimode when f r exceeds 3 GHz 
as obtained experimentally. Also, it is observed that tum-on jitter does not depend on off- 
states as long as the off-states is below threshold, but if the off state is above threshold 
jitter decreases. Additionally, it is found that jitter increases due to bit pattern effect when 
the data rate is high enough and pseudo-random pattern is applied. 



Chapter l 



The vertical cavity surface emitting laser (VCSEL) is emerging as the light source of choice 
for modem high speed, short wavelength communication systems.For short distance links it is 
of great interest to use Vertical Cavity Surface Emitting Lasers(VCSELs) as light sources[2] 
In order to simplify the driving circuits as well as to reduce the electrical power consumption, 
zero(or constant) biased operation is required. It is well known that this kind of operation 
leads to a significant turn-on jitter due to spontaneous emission and the bit pattern effect. In a 
digital transmission system an enhancement of the tum-on jitter increases the bit error 
rate(BER) especially at high data rate, which can affect the performance of the system 
significantly[2j. 


1.1 Vertical Cavity Surface Emitting Lasers (VCSELs) : 

An introduction 

The information below about VCSEL has been reproduced from the journal of Siemens Fibei 
Optics, Berlin, Germany [1]. 

Introduction 


The difference in the basic structure of a conventional edge emitting type laser and a VCSEL 
is that for the edge emitting laser the light propagation is in the plane of the wafer whereas it 
is perpendicular to the wafer surface in the VCSEL( Fig.l). A standard length for edge 
emitting lasers is 300pm whereas the active length in a VCSEL is typically only 24nm. This is 
the length where the lightwave must pick up all energy. Therefore mirrors with a very high 
reflectivity exceeding 99% are required in a VCSEL to enable laser operation,. In edge 
emitting laser the cleaved facets with a reflectivity of approximately 30% are sufficient to 
support lasing. Due to the waveguide design the output beam of conventional lasers exhibits 
an asymmetric shape with a relatively large divergence. In a VCSEL we have a large amount 
of freedom in the design of the shape of the active area. So we can tailor it to fit best to our 
application. In particular we can choose a round shape to obtain a highly symmetric circular 
output beam with a small divergence angle which eases coupling to optical fibers 

The good news for industrial applications is the potential low fabrication and packaging cost 
of a VCSEL. In fact, this is what drives the enormous interest in this young technology. The 
reason for low cost is that we can employ standard IC fabrication technologies to make the 
devices. All the fabrication and testing is done on wafer level with no cleaving and expensive 
handling of laser bars necessary. Since there is no cleaving, VCSEL wafers are not thinned to 
a fragile membrane like conventional lasers making the handling of these thinned wafers a 
nightmare for any production line. 



a Edge Em tt ng Lase 




Fig 1 


Basic difference between Edge-Emitting Laser and VCSEL 


History f 11 

The idea to built a semiconductor laser diode with a vertically oriented optical resonator wa i 
originally developed by Professor Kcnichi Iga from the Tokyo Institute of Technology in 
Japan in the late 70’s. In these lasers the light propagates perpendicular to the wafer surface 
within a very short cavity. The mirrors for optical feedback are made of dielectric multilayer 
stacks on either side of the epitaxially grown semiconductor active layer. After five years 
research Professor Iga’s group was able to demonstrate the first pulsed operating VCSEL at 
room temperature in 1 984. Four years later continuous wave operation at room temperature 
was achieved with threshold currents of 32mA. 

The major breakthrough in VCSEL technology happened one year later when Jack Jewell 
from AT & T bell Labs, New Jersey, USA, and Axel Scherer and Jim Harbinson from 
Bellcore, New Jersey, USA, together with some coworker -could demonstrate the first all- 
epitaxially grown VCSEL on a GaA-s substrate. Many hundred extremely thin semiconducto 
layers built the highly reflecting mirrors of this device. One quantum well with a thickness o 
just 20 atomic layers provided all the optical amplification of the travelling lightwave. These 
micro-cavity laser diodes showed sensationally low threshold currents of 1 .5 mA. 

These results initiated a boom in VCSEL research activities worldwide. Attracted by the 
remarkable research results companies began to show interest. Today there are eight 
companies and several research institutes in the U.S. with VCSEL activities. Some compann 
are manufacturing VCSEL already. The rest of the world follows at some slower speed but i: 
starting to turn the sails in the wind. 
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design must provide an optical confinement for the lightwave prop 
electrical confinement for the current which drives the laser. This i 
her laser. Among the various approaches to achieve these two main 
mostly favoured in fig2. The light in the VCSEL bounces up and d* 
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eedback from the mirrors leads to an optical intensity within the la: 
) orders of magnitude stronger than the electrical field intensity of t 
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laser current flows from a contact metallization on top of the etche 
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aplantation damage can affect the reliability of the laser. The optica 
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oneywell and a lot of other companies due to their relatively easy f: 


type employs the selective oxidation of a very thin AlAs layer in tl 
lighly resistive layer of A10 x just above the active layer. This oxide 
laser current and builds an optical waveguide at the same time Th 
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wavelengthf 11 

VCSEL work so far has been done on devices with an emission wavelength oi 
[i for data communication application. Table 1 gives an overview of the differs 
ed or discussed for VCSELs from the visible blue/green to the long wavelengt 
telecom applications. 

ate rial, performance and difficulties for VCSELs in the different wave len 
m visible blue/green to long wavelength infrared. fli 
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mirror technology, 
high temp.performance 


the wavelength range between 780nm and 980nm are based on GaAs or InGa/ 
11s and are relatively easy to fabricate. They exhibit excellent performance and 
n at various places already. Visible VCSELs with an emission wavelength of 
i have been made using InAlGaP active quantum wells. They will come into 
oon. The main application are data storage, plastic optical fibers and laser 
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From Table 2 vve see that VCSELs combine the best of LEDs and laser diodes They are as 
cheap and as easy to handle and to package as LEDs but they provide the speed and the high 
output power of a laser diode. In fact it is HP s philosophy to treat VCSELs simply as fast 
LEDs. It is obvious that such a device is very attractive for many application. 


Advantages of VCSELs 

Performance issues 


ffl Surface-Normal output 
H Ultra-low lasing threshold 
S Dynamic single-mode operation 
H Circular, low divergence astigmatic-free beam 

■ Extremely high electrical to optical power conversion efficiency 

■ Thermally stable operation 
M High speed 

M Extremely small active volumes 
B 1- and 2-D arrays 

Systems issues 

H Efficient coupling to fibers 
1! Insensitivity to feedback 
H Densely packed two-dimensional arrays 

■ Long lifetime 

BB Compatibility with electronics 
Manufacturing issues 

Bl Low cost wafer-scale chip production 

■ On-wafer testing and screening 

■ Simplified mounting and packaging 
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applications of vertical-cavity lasers and laser arrays. 


few examples of the application for VCSELs. In serial optical links VCSELs 
> due to the higher speed. Gigabit Ethernet is one example for this kind of system 
are used in parallel optical data links aim at the server, supercomputer, and 
;h market where very high data rates have to be transmitted over distances up to 
gher speed of VCSELs compared to LEDs is also attractive for high speed opto 
i bit rate of more than lOMbits/s. The volume for opto couplers is high comparec 
nn market, but the price per VCSEL must be very low in order to compete on the 
:a3 free space communication systems work very much like a TV remote control 
rovides data transfer for example from a PC to a printer through air by IR 
ead of using a cable. 

are a huge market for VCSELs. In the traditional way of a CD player the 
s encoded optically on the disk and the laser beam reads the information. 

high power capability and high efficiency of VCSEL arrays is attractive for 
d state lasers. Furthermore it is a lot easier to attach a VCSEL than an edge 
r to the solid state laser. This reduces the cost of the whole system dramatically 
larket is full of different application. Large market volumes are predicted for 
mg in automotive applications and encoders to find the exact position of rotating 
parts. Also laser printers which mostly LEDs today are a possible market for 
, arrays. 

ications VCSELs are just a replacement for conventional lasers or LEDs like in 
iks where VCSELs are the better choice due to lower cost and/or higher 
Other application require some of the unique features of VCSELs and it is very 
e job done with conventional lasers. This is true for most applications which 
arrays like disk sensors, high power pump sources, and laser printers. 
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(_ onclusiou 


VCSELs arc very attractive for the optoelectronics industry due to their specific feature! 

1 . Potentially low fabrication and packaging cost ( non hermetically sealed packages), 

2. VCSELs are attractive sources for Sow-cost, high speed and high-performance, 
manufacturable parallel optical interconnect modules. 

3. Standard 1C fabrication processes on full wafers ( no cleaving, no handling of small 
pieces). 

4. Wafer scale testing. 

5. Single mode emission. 

6 . Circular output beam shape with low divergence. 

7. No HR/AR facet coating. 

8. Easy fiber coupling, no additional lenses necessaiy. 

9. Easy integration with other optical elements. 

10. VCSELs can also be fabricated with different emission wavelengths within a 2D 
array in a controllable fashion. Thus, wavelength division multiplexing (WDM) 
applications arc feasible. 
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Introduction 


It is well known that biasing a laser near or below threshold can significantly increase the bit- 
error rate (BER) of optical communication links operating at gigabits rates[2]. However, there 
are application such as optical interconnects, where zero biased modulation is preferred since 
it eliminates the need for optical monitoring and feedback control of the bias point. It is of 
great interest to operate VCSELs in optical interconnects with zero bias since the laser driving 
circuit is simplified and power consumption can be reduced[2]. 

Considering the corresponding tum-on event, the BER degradation due to an improper biasing 
can be easily understood. When a laser is suddenly switched from below to above threshold, 
the optical emission will be randomly delayed. The turn-on delay is defined as the time for the 
photon population to grow from its initial value at the moment of excitation to the stationary 
value of the on state [2]. Due to spontaneous emission during tum-on there is random 
fluctuation of tum-on delay. The turn-on jitter is defined as the standard deviation of the 
probability density function (PDF) of this turn-on delay [2]. There are two main mechanisms 
which contributes to the tum-on jitter: spontaneous emission and bit-pattem effect. Due to this 
turn-on jitter there is significant BER at gigabit data rates. Although at low data rates the 
BLR is not significant due to this jitter[2j. 


Theory 


An analytical expression to describe the probability density function (PDF) for the tum-on 
delay for a single-mode edge emitting laser was given in [2]. In contrast to the single TE- 
polarization state in edge emitting lasers usually two orthogonal polarizations contribute to the 
emitted light of single mode VCSELs. The PDF for the tum-on delay of a transversal single 
mode VCSEL with two orthogonal polarizations can be derived in a similar way as shown in 
[3], 

Let us consider a tum-on event as sketched in Fig. 2.2. 1. The laser drive current is switched on 
at t=0 from its bias value d 0 ff<I th to l on > hh w dh hh denoting the threshold current. The 
photon number S within the cavity will then increase starting from the initial value S(0). In 
order to calculate the statistical properties of the tum-on delay t d a model has been used 
here(e.g., [2], [3]) which splits the behaviour of the laser into a stochastic and a deterministic 
regime similar to that originally proposed in [3J. In contrast to this approach, it has been 
assumed here that for every tum-on event (he., for every combination I 0 ff, ion) there is a 
certain t c which separates the two regimes. For t<-t c the photon number is very small and 
fluctuates strongly due to the dominance of the spontaneous emission.. This leads to an 
exponential PDF of the photon population which can be assumed is still valid at the crossing 
time t c according to ^ ^ 


P s ) (ys„) C x P ( s/s ( 22 i)«»*kA j.280 


Photon Number 



where, S c , (S c ) denotes the photon number and its average at t-t c , respectively. (S c ) 
corresponds to the ^absorbing barrier" as introduced in [4], [5]. In contrast to the crossing 
time t c , the value of (S c )is independent of the specific turn-on event. For t>t c , the laser enters 
the deterministic region, where the additional time t cn =t d -t c for the photon population to 1 
grow from S c to the stationary value at the on-state S on is given as [6] 


= (]/2irf, )V2In(S m /S t ) ...(2.2.2) 


with/ r denoting the relaxation resonance frequency. 

Since VCSELs have two polarizations, the PDF of the photon number in each polarization 
i— 1,2 of the SM VCSEL is given by 

PfS'j) = (l/{s cJ ))exp(-S tJ /{S CJ )) , i=l,2 ...(2.2.3) 


At the crossing time r c , S t -,i and (s cj ) denote the photon number and its mean value in the 
considered polarizations. 

In the case of the SM-VCSEL both polarizations contribute to the total photon number 
S c = S cA + S c2 and its mean value {S c ) = (s cd ) + {s c<2 ) . 

The PDF of S c is given by the convolution of the PDFs of both polarizations yielding 


P(SJ = (*Sj{S'Y) • exp(-2 SJ(S,}) ...(2.2.4) 



where (5 /, ) - (S c>2 J = (S c )/ 2 is assumed. 

In the deterministic region t on is given by eqn, (2.2.2) for a known S c . 

The PDF .or the photon number 5,. can be transformed into a PDF for the tum-on deiay usin* 
the relation p(t on )dt - p(S c )dS c yielding 


p(t °"> 4co '- 2t '>^ S o*'7( S c} 2 )- ex PE'-(® r t on ) 2 ]- exp[- 2S on /(S (; }exp[- (ca r t on )72]3 ...(2.2.5) 

where co r = 2nf r . 

The turn-on jitter is given by the standard deviation of eqn. (2.2.5). 

Once we have determined the value of the “ absorbing barrier 4 ' (Sj (which varies from one 
laser to other), the probability distribution according to (2.2.5)w only dependent onf. and S on . 
It turns out that p(t on ) depends weakly on the photon ratio S ort /(s c ) compared to the 

dependence onf r . Furthermore, since/, is proportional to , the relaxation resonance 

frequency can be considered as the main relevant parameter. 



t - < t > f [ ps ] 


n 



00 s 



Fig. 2.2.3 



Fig. (2.2.2) and (2.2.3) shows the PDF of the turn-on delay for f r =3.7 GHz and/ r =2.985 GHz 
respectively with different (S c )s using the eqn. (2.2.5). The turn-on jitter for different (S c )s 
versus relaxation resonance frequency f r is plotted in fig. (2.2.4) using the eqn. (2.2.5). So, 
according to the theoretical model above, for a particular VCSEL(i.e for a given (s,) ) the 
jitter decreases with the increase of f r i.e, with the increase of on-state level. 
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1 3 General Normalized Rate Equations 


Starting from the basic rate equations the normalized rate equations has been derived here 
which is suitable for solving it numerically [7]. 

We have the basic rate equation for photons as 


dS/dt = S(R„-l/t ph ) + R Ip ...(2.3.0) 

Where spontaneous emission coefficient /? =r> and stimulated emission coefficient R st = 
Tr Jf and S is the photon numbers. 
r sp - Spontaneous emission rate. 
r tl ~ Stimulated emission rate. 

F *= Confinement factor. 
t ph= Photon life time . 

In the rate equation above the term SR s , defines the number of photons generated per unit time 
due to stimulated emission while the term S/t P h defines the number of photons lost per unit 
time due to recombination. The term R sp defines the spontaneous generation of photons. 

For laser operation we can introduce a factor called inversion factor n ip . Under typical lasing 
condition n rn - 25 

sp 

After introducing n sp we can write 

Rsp= n sp Rs'^Hsp^pk .-(2.3.1) 

The expression in eqn. (2.3.1) is valid for I>I(b.Thus we can neglect the photon number 
produced by spontaneous emission and we get for stimulated emission R st (n = n th ) => 1 ft ph 

Now we have the rate equation for.electrons as : 


dn/dt - (I/e) - R(n) - R J( S ...(2.3.2) 


where, the spontaneous recombination rate is represented by R(n ) , the injection current by / 
and e is the elementary charge. SR st defines the carrier consumption due to stimulated 
emission. 

For the recombination rate R(n ) we can assume R(n)~n/ t r .where n is the carrier number 
in the laser cavity and i f is the carrier life time in the active region. 



Now for a particular photon energy gain depends on the earner density as 


g st (n) = (a/V) (n - n Q ) ...{2.3.3) 

Where g st ( n) is the stimulated gain at the carrier density n IV and n Q IV is the earner 
density at population inversion , a is some constant and V is the volume of the active region. 

Taking the confinement factor into account we can write effective gain g = Vg s! . 

Now , since stimulated gain coefficient g n = r st jv $r , g = R st /v gr 
Where, v gr is the group velocity of the optical wave in the active region. 

So, the threshold gain g th is given by 

g(nj = R st {n th )jv sr =1 /t ph v sr = g Ih ...(2.3.4) 

So, we can write normalised gain G= g/g lh = R st t ph ...(2.3.5) 

Since, n ~ n lh above threshold, we have g ~ g lh i.e, G*1 for lasing operation. 

Now, the condition for stationary laser operation is 

~a s )L) = 1 ...(2.3.6) 

where, a s is scattering loss in the cavity per unit length, L is the length of the cavity and Ri, 
R 2 are the reflectivities of the two mirrors. 

Therefore we can write 

g th =«,+ (1/2L) ln(l//?iR 2 ) ...(2.3.7) 

Now, n th is the threshold carrier number, where the loss in the cavity is compensated by the 
stimulated emission and is given by 

n lh =n 0 + (l/(t ph v gr T(^V))) ...(2.3.8) 

By introducing the normalized gain G= g / g lh = R st t ph we can write the rate equations 
(2.3.0) and (2.3.2) as follows 


dn/ dt — I/e — R(n) — GS j t ph ..,(2.3.9) 

dS/dt = (S/t ph ) (G - 1) + K , ot R sp ...(2.3.10) 


where K is the enhancement factor for spontaneous emission. 

l£> 


In order to account for the nonlinear gain the gain G can be splitted into linear and nonlinear 
part as 


G = G i (1-k p P) = G i (1-k,S) .42.3.11) 


where G, is the linear gain and 

K is the gain compression coefficient related to optical power and 
K s is the gain compression coefficient related to photon numbers. 

The relation between k s and k is given by 

K s = {\f2)hw gt - (In (l/i?,i? 2 )) -K p /(2L) ...(2.3.12) 

where h = Planck’s constant 

v = Freq. of the emitted light 

In order to solve the above differential equations numerically it is necessary to normalized the 
equations as 

S N =S/(n th t ph /t e ) 

N = n/n th ...(2.3.13) 

In = VK 


where S N , N and f N are normalized photon number, normalized carrier number and 
normalized current respectively and I lh = n, h ejt e ...(2.3.14) 

Now the linear gain G t (n) » G l {n th ) + {dGjdn)(n--n th ) ~ 1 +(dG l /dn)(n — n th ) ...(2.3.15) 

So, from (2.3.13) and (2.3.15) the linear gain G, can be written as 

G t ~l + (dGjdN) (N - 1) ...(2.3.16) 

where, dG,/dN = (n lh f g th )T(ajV) ...(2.3.17) 

After introducing eqns. (2.3.13), (2.3.14) in eqns. (2.3.9) and (2.3.10) finally we get the 
normalized and dimensionless rate eqns. for electrons and photons which are given by 


dS N /dt= (l/t ph ) ((G - 1)S N +P S ) 
dN/dt - (l/t e )(I N -N - GS n ) 


.(2.3.18) 


where = K lot n sp tJ(n [k t ph ) 
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'Hie gam can be written as 


G (N ) = [1+ (dGj/dN)(N - 1)][1 - kS n J ...(2.3.19) 


where, 

= Normalized spontaneous emission rate. 

K — Normalized gain saturation coefficient. 

k is related to k s as 

K = K,n t h t ph jt e ...(2.3.20) 


Simulation of noise sources : 


The noise due to spontaneous emission is realized by the Langevin noise sources F n and F s 
given by [8] 


F s (f, ) = Q2Sy (« w )<„ ft /VAf ) • X, . 

F„ (/, ) = (ft /Vv?) *. - ) ■ x. 


where At is the distance between the discreet time slots t, and t,_/. 

The random variables x e and x n are gaussian distributed as 

<^>=<o=° 

<o=<o=l 

The random variables are not correlated to each other. 

After inserting the noise sources in the rate equations we can obtain final rate eqns. as 


dS N fdt = (1 jt ph ){(G - 1)S N + A ) + F s (f, )/t ph 
dN/dt = (ljt e )(I n -N-GS n ) + F n (f { )jt c 
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If the power output from the 1 st mirror(output mirror) is PSP I and the power from the 2 
mirror is PSP2 then 

the total power Ptotal - PSP1+PSP2 

where, PSP2 = ((I - i? 2 )/(l -RJ)- ^RjR^-PSPl ...(2.4.1) 

So, the number of photons generated in the cavity corresponding to this power Ptotal is 
SP = (2 ■ Ptotal • L • l)/(h v gt c- ln(l jR x R 2 ))) ...(2.4.2) 

where, A is the wavelength of emitted light 
h is Planck’s constant 
c is velocity of light 

The normalized photon number is 

SN =%,,<WoH 2A3) 

where n lH *( t ph jt e ) describes the photon number at 1=21 

and the normalized current for this photon number can be written as 

INbias = (1 + SN) ...(2.4.4) 

So the bias/on current is 

Jbias =I th (l + SN) ...(2.4.5) 

The bias current above is the on state current when the modulation current is applied. 


LaseitVCShL) 


i 


For VCSEL the light propagation is in the perpendicular direction to the wafer surface, The 
light in the VCSEL bounces up and down between two highly reflecting multilayer Bragg 
reflectors on either side of the active layer. Most of the VCSELs are round shape to obtain 
highly symmetric circular output beam. 

So, wc have for VCSEL the radius =r 

The effective length = L e ff 

The thickness of the active region - d 

So, the volume of the active region V — n.r 2 d ...(2.5.0) 

Here we have considered the single mode single polarization VCSEL without the effect of 
multilayer system. 

For VCSEL we have two confinements in longitudinal and transverse directions. In transverse 
direction the confinement factor is taken as F and in the longitudinal direction the 
confinement factor is taken as F f . Since the thickness of the active region of the VCSEL is 
very small compared to edge emitting laser we have r z very small for VCSEL and because of 
full confinement in the transverse direction F = l. 


For VCSEL we have d/L eJ} * ...(2.5.1) 


So, for VCSEL the effective gain 
...(2.5.2) 

g = rr t s„=rr>/iO(/.-n.) 

[ from eqn. (2.3.3)] 

The condition of stationary laser operation is 


...(2.5.3) 

exp((g-a,)L^) = l 

[from eqn. (2.3.6)] 

So that for VCSEL we can write the threshold gain as 


...(2.5.4) 

go, =<*,+( l /2L, # ) in( !/ /? 2 ) = 

= rr I (a/v)(«.-B 0 ) 

and 

n, h =n 0 +{\j(t rh vJY z {ajV))) 

...(2.5.5) 

Also, 

dG, /dN = (n sh /g,jiT\ (ajV) 

...(2.5.6) 

and 

k s ={\j2)hw gr - (In(I /R t R z ))-K 

„/(2L tf ) ...(2.5.7) 
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For VCSEL the photon number in the cavity corresponding to power Ptotal is 


SP = (2 * Ptotal • L tff • A) /(h • V gr • c • Sn(l/ R \ R 2 ))) 

These are the changes of equations which has been incorporated in simulation program for 
VCSEL. 
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the Simulation of Turn-On itter 


2.1 Introduction 


We have at our disposal two rate equations with noise sources and the expression for non- 
linear gain as follows. 


G(N) = [l + dG/dN(N-l)][l- kS n ] ...(3.1.0) 

dS„ /dt = (^t pt )((G - 1)S N + p, ) + (F s (t,)/ t5h ) 

dN/dt = (l/t t ) (I„ - N - GS „ ) + (F„ (t, )/t, ) "' ( ‘ ' 


The numerical solution of the coupled differential equations given by (3.1.0) has been done b) 
4 h order Runge Kutta method without adaptive stepsize control using the software 
Lab VIEW. 

The random numbers x e and x n in the Langevin noise sources F s and F n are generated by 
standard gaussian noise sources from Lab VIEW standard tool. These Langevin noise sources 
are added to the rate equations to simulate the spontaneous emission noise. 

A subprogram which generates the necessary parameters for solving the rate eqns. has been 
developed. This parameter file needs some basic parameters of the laser as input to give some 
output parameters which are taken as the input of the main program. 

The laser characteristics graph has been plotted first by changing the power output and 
observing the corresponding bias current. 

Initially the bias current is applied to make the laser in steady state after few oscillation. This 
bias current is the on-current I on when modulation current is applied. For a particulanoff - 
current I off different on-current i.e different output power P on at the on state of laser is applied 
Also for a particular on-current different I„ (f is applied to observe tiie effect of off-current 
First the 1010 pattern is applied and then the pseudo random pattern is applied to observe the 
bit pattern effect. 

For simulation with the VCSEL provided by Prof. Ebeling, University of Ulm, Germany, the 
data rate taken for simulation of jitter is 160 Mbits/s to match with the experimental results. 
Also to observe the effect of high data rate on jitter the simulation is carried on for / Gbits/s. 
The duration of a time slot ( i.e step size for numerical solution) Af =le-12 is taken for 
convenience. The total number of bits taken for simulation are 5000. The photon turn on dela 
is ealeu ated whenever the current changes from off state to on state So during the time whei 


5000 bits are transmitted many turn on delays (t on ) are obtained The probability density 
function (PDF) of these t on s is plotted The mean t on is calculated and the RMS of this PDF 
gives the simulated jitter. 

Also the average relaxation oscillation frequency f r is measured automatically in the main 
program by using various Lab VIEW tools. 

In simulation the values of jitter and f r are calculated for several on states i.e for several 
output power P on with zero bias and compared with the experimental results as well as with 
the theoretical model[ 2] and then the simulation is carried out with different off-states. 

Then the graph between rms tum-on jitter and relaxation resonance frequency is plotted to 
compare with the theoretical model as well as with the experimental measurements. 

In this simulation of turn on jitter and calculation of f r for Ebeling VCSEL some parameters 
have been chosen logically apart from some standard supplied parameters of the laser to 
match with the experimental results. 

The parameters which are supplied for Ebeling VCSEL are radius( r), length of the active 
region ( d ), effective length (L e jf ), wavelength of the emitted light ( X ), group velocity( v gr ). 

The parameters which are chosen are electron life time(t e ) , photon life time ( t ph ), reflectivity 
of the mirrors ( Ri,R 2 ) so that a good fit with the experimental results and theoretical model 
can be achieved. 


This type of simulation is called simulation verification which approximately represent the 
real world system. In this simulation model some simplifying assumption has been made e.g, 
single mode single polarization VCSEL without considering the multiple quantum well effect. 
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escn 


diagram of Che program 


tgram of the program which has been developed by LabView is shown belc 


Input 

Parameters 


Modulation current 


current 


Parameter subprogram 
(Calculation of output 
parameters and initial 
values For numerical solution) 


Total current 
♦ 1010 

$ Pseudo random bits 


Initial values of 
photon number & 
electron number 


Gaussian white 


(noise 

subprogram) 


K-K sub program 

* Normalised rate equations 

* Equation for gain 

Runga Kutta method - one step 
(Numerical solution of the rate eqns.) 


ace initial values 
lext values 


Next values Of 
photon number and 
electron number with 
time step At 


Plotting of 
output waveform 


Simulation 
complete ? 


Calculation of f and 
RMS jitter. 
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2 2.2 The parameter subprogram 


When the main program is runned, it calls the parameter subprogramfpararaeter. 
that user can give necessary input parameters of the particular VCSEL to calculate s< 
derived parameters which are necessary for solving the rate equations numerically ai 
simulation of jitter. The input parameters of the Ebeling VCSEL which have been u: 
LabVIEW subprogram parameter.vi for the simulation are given below: 

M Output power through mirror I (PSP 1 ~P on) -variable 
B Wa velength(lanula)—870nm 
& Radius of the VCSEL(rad)=2pm 
M Thickness of the active region(d)=24nm 
M Effective length(Leff)~3.5prn 

& Group velocity of the optical wave in the active region(vgr)~8.3e7m/s 
B Photon life time(tph)~4.4ps 
B Electron life time(te}~ 1.2ns 

B Carrier number at population inversion(n0)=6.6e5 
Mf Reflectivity of mirror }(R1)=0.99 
MS Reflectivity of mirror2(R2)=0. 99 
B Gain proportional coefficient dgstdn=a )=2. Qe-20tn 
M f Time slot interval Step size)(dx)~Ips 

SB Gain compression coefficient related to optical power(kp)~2SW J 
SB Longitudinal confinement factor( gamma- T z )=6.8e-3 

With these input parameters some output parameters are calculated which are neces 
the main program to run. These parameters are given below: 

M Normalized bias/on current corresponding to the power PSP I ( INb ias) = variable 
B Normalized photon number in on state corresponding to the power PSPI(SN-on, 
value of normalized photon number needed for simulation(SN)— variable 
SB Normalized value of initial electron number(N)-l 
SB Photon life time(tph)~4.4ps 
SB Electron life time(te)- 1.2ns 

SB Starting value of independent variable! time) for numerical solution(x)~0 
SB Time slot interval Step size)( dx)=Ips 
SB Normalized gain saturation coefficient(k)-1.666e-2 
B Normalized spontaneous emission rate(betas~ (i s )=8.9933e-5 

K dG/dN- 1.1087 

B Threshold carrier number(nth) =6. 7323 e6 
B Threshold gain(gth)~2.7382e3m 1 

The above derived parameters are calculated in the parameter subprogram accordm 
equations described in chapter 2.3. 

When the main program is executed first the parameter subprogram runs in a numb 
sequence to calculate output parameters which are subsequently fed to the R-K 
subprogram(newRK4.vi). Initially the laser is assumed to be in steady state with t 
powcr(PSPl) so that the initial (steady state) value of electron number (N=l) and p 
number(S N-on) needed tor numerical solution ot the rale equations me calculated a 
the R K subprogram So the bias to the aser is given from the parameter file 
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The VCSEL supplied by Prof Ebelmg has multilayer system For emission near 870nm 
wavelength the one wavelength thick inner cavity contains three each 8 nm thick GaAs 
quantum wells. So the thickness of the active medium is taken as d=24 nm. Since longitudinal 
confinement factor F 7 ~ dj L eff and 2^=3. 5pm, so T z = 6.8e-3 has been taken correctly. The 

radius of the VCSEL is supplied as 2 pm. In the simulation model of the VCSEL the 
experimentally obtained threshold current I t h =0.9 mA and slope of the laser characteristics 
( dP/dI=0.775 ) above threshold are matched first by the following formulas 


fh ~ n t h e{t e ...(3.2.2. 1) 

where n th = n 0 + (l/(t p v gr T z (ajV))) ...(3.2.2.2) 

and dP/dI = (hv gr cln(l/R l R 1 )/4L eff te)-t p ...(3.2.2.3) 

where R 1 = R 2 = 0.99 and t p = 4 ps , t e = 1.2 ns are taken. 

A separate program (fr.vi) has been developed to calculate Ibias, f r and l th by the following 
formulas: 

Ibias = ( n 0 e/t e ) + ( e /y e v gr F J ( C i/V)) +(4e(FSPl)L # A/f/v gr cln(l/i? 1 ^)) ,..(3.2.2.4) 

= C l /t e + C 2 /t p t e+ C 3 (PSPl)/t p 
[from equations (2.4,1) to (2.4.5) in chapter 2] 

C, = n Q e — 1.056c — 13 

where C 2 = e(y g X t («A0 = 4.29 e - 24 

C 3 = 4eL eg Xlhv ?r c\n{\jR,R 2 ) = 5.9c - 12 

and f r ~ (1/2 7t)-^jT z v gr (a/V)/e Ibias- ( n tk ejt~ ) ...(3.2.2.5) 

= C 5 [^] Jbias -n lH ejt t ] 

[ from the basic equation f r = (ifltz) ■ ^(1 fn th ){dG/dN) (I bias - 1 th )fet ph ] 

where C 5 = (l/2n) ■ J T z v gr (ajV)/e = 7.68dO 
and I th is given by the formula (3.2.2. 1) using n tfl (3 .2.2.2). 

In the program fr.vi a range of values of t P h (from 3.5 ps to 4.4 ps) and a range of values of t e 
(from 0.5 ns to 2.3 ns ) has been used to calculate f r Jbias and I th for a particular power PSP1 
and those values of t ph and t e are taken finally which gives a good fit with the experimental 
values of f r Jbias and I t h- 

Finally the values of photon life time t p h and electron life time t e are taken as 4.4 ps and 1.2 ns 
respectively with which the simulation results are matching well with experimental results. 

K p is chosen as 25 for sufficient damping to limit the oscillation. 
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12.2.3* The noise subprogram 


In the noise subprogram (noi.vi) the Langevin noise sources F n and F s are generated to add 
these functions to the rate equations in the main program. The random variables x e and x n are 
generated in the main program by standard LabVIEW gaussian noise generators with standard 
deviation equal to 1 and these are fed in the noise subprogram. For every step of solving the 
rale equations numerically one value of noise is generated and applied to the rule equations. 
These noise sources are added to simulate the spontaneous emission noise so that one can 
observe fluctuation of turn on delay and jitter can be simulated. 


; 2.2.4 The R-K subprogram 


The R-K subprogram(newRK4.vi) consists of the following subprograms : 

■ newfxy.vi 

■ newarr.vi 

H newdSNdt_dNdt.vi 
a newGN.vi 

In the R-K subprogram the 4 h order Runge-Kutta method is applied to solve the rate 
equations simultaneously. In this subprogram one step R-K method has been done to calculate 
the new values of photon number(SN) and electron number(N) after giving the initial values 
of these using step size Ax = Ips . The R-K method can be described as following: 

If we have a set of N first order ordinary differential equations as[9] 

dy,(x)/dx = Mx, y lt y N ) where 1=7,2 N 

Usually, it is the nature of the boundary conditions that determines which numerical methods 
will be feasible to solve the equations. In the initial value problems all the y. are given at some 
starting value and it is desired to find the >*,•’ s at some final point x f , or at some discrete list 
of points. 

One practical numerical method for solving initial' value problems is Runge-Kutta method 
according to which for any equation if x n and y n are the previous values and if the value of x R 
is increased by step size h U, 1+1 = x„ +h ) then we can calculate the new value of y n as 

y fl+ , -y, +KJ6+ K 2 j3+ KJ3+ KJ6 
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where, 


^ h f(x m y ) 

A h f(x + h/2 v +K /?) 
^3 ~ h‘ f(*„ + h/2, y n 4- A't /2) 
^4 =h- f (x n + h, y n + K z ) 


1 he 4 order Runge-Kutta method requires four evaluations of the above right-hand side per 
step A. 

Here only two rate equations has been solved, so N~2 and yt~SN and yh~N , where SN and N 
are normalized photon number and normalized electron number respectively and h — Ax . 

In the subprogram newGN.vi the equation for nonlinear gain G(N) (eqn.2.3.19) is constructed 
by standard LabVIEW tools which is required in the rate equations. 

In the subprogram newdSNdt_dNdt.vi the rate equations for photon and electron has been 
constructed with the noise sources added from the noise subprogram and the bias current as 
well as the modulation current is applied here from the main program. The newGN.vi is a 
subprogram of this program. 

In the subprogram newfxy.vi the photon number SN and electron number N have been put in 

an array and applied in the subprogram newarr.vi which separates the elements in the array 

and the separated elements (SN,N) are then put in the subprogram newdSNdt_dNdt,vi. 

* 

In the program newRK4.vi the program newfxy.vi has been used as a subprogram to solve the 
rate equations numerically by R-K method. 

2-2.5; The main program 


The main program VCSELPRS2.vi consists of the following subprograms. 
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When the mam program is executed it takes first the parameters necessary from the paramete 
file to solve the rate equations Initially the laser is biased with some bias current 
corresponding to the power PSP1 and after the laser comes to the steady state the modulation 
current is applied from the main program where the on-state is taken as the bias level. The 
simulation is carried on for several off state. Two types of modulation current is applied -1010 
pattern and pseudo-random pattern, The subprogram PRB.vi generates the pseudo-random 
bits taken from the data file PRB.dat. From the front panel of the program the data rate is 
given and corresponding bit duration is calculated as well as the number of sampling points 
during one bit period are calculated. The program is executed for several thousand bits and 
the variation of normalized carrier numbers and normalized photon numbers are plotted with 
time with the modulation current. In every sampling step of numerical solution of the 
equations one value of noise from the gaussian noise subprogram is added to the equations. 
Also the turn-on delay is calculated when the laser is switched from off-state to on-state and 
the turn on delay is taken as the time when the photon number reaches the steady state value 
first during tum-on event. The mean turn-on delay as well as probability density 
frmction(PDF) of the turn-on delay is calculated in the subprogram PDF-Direct.vi. The rms 
jitter is also calculated in the same subprogram. The relaxation resonance frequency f r is also 
calculated in every on-state and the average f r is calculated. The main program is executed for 
several on-state power PSP1. 
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T 3 


ition results and conclusions 


haracteristics 


ymg the modulation current the characteristics of Ebeling VCSEL is it 
g the parameter file. After introducing the parameters of the VCSEL v 
mtput power vs. bias current by changing different power PSP! [Table 
ics graph obtained is matching well with the experimental graph as sht 
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Current(mA) --> 


4.1.1 


P-I Characteristics of Ebeling VCSE3L 



4.1.2 


Experimentally obtained P-I characteristics of Ebeling VCSEL. 






mu it i on it s assumed t nt the VCSEL has on y one po anzation But from th 
nt s oil isuiunui s (fig 4 l 2) it h is been foun 1 1 1 1 ic Fhc mg VCSFI h is 
sons- one dominant and other is suppressed polarization which has significant 
mg turn-on time under zero-bias conditionfor below threshold)! 10]. But from 
i islics graph (tig 4. ! .2) obtained experimentally it shows the power of the sup 
ion is very small compared to the dominant one. So although in simulation it i 
that there is only one polarization , the simulated P-I curve may be matched v 
experimental one. 


[inflation of turn-011 jitter 

to investigate the turn-on behaviour of the VCSEL the modulation current is a 
e experimental measurements have been done with bitrate 160Mbits/s with z 
I 1010 bit pattern, the same data rate has been used in simulation also to com 
measurement results. Additionally, the simulation has been done for differeni 
observe the effects of off-state on jitter. The sampling time interval for solvin 
lal equation is taken as 1 ps. 


Table 2: Results obtained by simulation for zero biased VCSEL 


Data rate=l 60 Mbit/s 
Bit pattern-1010 


m 

fr(GHz) 

measured 

fr(GHz) 

simulated 

JlTTER(ps) 

measured 

JITTER(ps) 

simulated 

7.5OE-04 

2 

2.18 

25.39 

22.34 

8.0QE-04 


2.25 


21.3 



2.44 


19.69 

1.1 IE-03 

2.77 

2,64 

17.16 

17.16 

1.34E-03 

2.94 

■Hun 

13.19 

15.5 

1 .50E-03 


3.06 


14.4 

1.88E-03 

3.33 

3.48 

7.78 

12.6 

2.29E-03 

3.7 

3.96 

5.75 

11 


shows the measured rms turn-on jitter vs. the relaxation resonance frequency 
with the results found by numerical simulation of the rate equations for zero l 
1010 pattern) at 160 Mbit/s [ Table 2]. It can be seen, that the turn-on jitter dec 
asmg values of/ r . By comparison with eqn, (2.2.5) in chapter 2 , we find that i 

g barrier" (S c ) for this specific laser is about 188. As evident the simulation i 

iod agreement with the measurement (/)■< 3 GHz) and theoretical model. Abov 
at sufficient high on-state the laser is multimode as obtained experimentally a 
tion of jitter vs / as obtained from measurements is not match ng with the s i 















e results from the simulation is matching with the theoretical graph smc 
Lodel is derived only for single mode VCSEL. Above f r =3 GHz the exis 
r ersal modes yields a reduction of tum-on jitter due to the increased degi 
the VCSEL to reach the stationary value at the on state. However, if the 
, the mode partition noise must be taken into account yielding an increa, 


i Jitter vs. fr for VCSEL (Ebeling), 4pm, 870nm 
160 Mb/s, Bitsequence: 1010 
\ 



1 1 f . 1 > 1 • 1 > 1 > j 
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Relaxation Resonance Frequency [GHz] 


Fig 4.2.1 


Turn-on jitter vs. relaxation resonance frequency 


, the dependence of turn-on jitter on the off-state has been investigated 
delay pdf according to eqn. (2.2.5) should be independent of the off-state 
This has been confirmed experimentally for edge-emitting laser[2]and f 
umulation result shows the same behaviour as it is expected. The jitter i ( 
as the off-state is below threshold. 

the turn-on jitter decreases as obtained from simulation results as it is ex 
ulated emission is dominant above threshold compared to spontaneous 
3le3]. 


Table 3: Simulation results for different off-states 



JITTER(ps) 

simulated 

INOFF=0 

JlTTER(ps) 
simulated j 

!NOFF=0.9 

JlTTER(ps) 

simulated 

INOFF=1.2 

7.50E-04 

22.34 

22.41 j 

12.5 

8.00E-04 

21.3 

20.64 

11.5 

9.50E-04 

19.69 

19.06 

10 

1.1 IE-03 

17.16 

16.74 

8.6 

1 34F-0.3 

15 5 1 

16 63 

7 3 
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i resu, Ls s...,ws tL.t relaxation resonance frequency Ji.es nu depends oi 


Cal culation of relaxation resonance frequency for different off-st ate 


rnM j 

fr(GHZ) 

simulated 

SNOFF=0 

fr(GHz) 

simulated 

!NOFF=0.9 

fr(GHz) 

simulated 

INOFF-1.2 

.50E-04 

2.18 

2.18 

2.4 


2.25 

2.25 

2.48 

.50E-04 

2.44 

2.44 

2.69 

.1 IE-03 

2.64 

2.64 

2.89 

.34E-03 

2.91 

2.9 

3.17 


3 

results are in good agreement with the measurements and with theorem 


GHz the single mode simulation results is not matching with the 
al results since the VCSEL is multimode. 

, independent of off-state as long as the off-state is below threshold. 

ate is above threshold the jitter decreases as expected. 

resonance frequency does not depend on off-states 


of bit-rate and bit pattern on jitter 


uilation of iitter at different data rates and bit patterns with zero bi 

tlM-.iiiini.iimr— , I T * I -V— IT r\ / . \ 


JITTER(ps) 
160 Mbit/s 
1010 pattern 

JlTTER(ps) 

160bit/s 

pseudo random 

JlTTER(ps) 

1 Gbit/s 

1010 pattern 

JlTTER(ps) 

1 Gbit/s 
Pseudo-rand 

27.52 

27.4 

26.12 

4i 

21.3 

21.81 

21.97 

91 

17.16 

17.17 

18.42 


14.4 

14.35 

14,52 


11 

11.17 

11.11 
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VCSEL(Ebelmg) 
Data rate = 160Mbits/s 
lamda=870nm 
1010 pattern 
loff=0 












MmiK 0 Pa " em i ‘° 0bSOTe * e behaviour of jitter on data rat, 
eotmng almost same as that of when the data rate is 160 Mbit/s[Table 

a PP Iied * 1 GWts/s due to bit pattern effect the jitter 

. 4^f \ 5 t u- PDF SrapK 0btained U not thc same nature as 

fora . ™ m 5 15 ^ ecause ,he mm -° n de,a y on the numbe, 

fore a one(l) comes. For pseudo random bit pattern this happens 

e°u r "m Ur a ° f ? PDF graph ° btaincd (F « «- 2 > » quite different 
from the dtstnbution graph that the probability of a single zero bef 

i the probability so that two zeros come and so on. The values of jittei 
i ion resu ts o not give the actual value of jitter since the PDF graph h 
i gives the measure of fluctuation of tum-on delay. So qualitatively w< 
jitter increases for pseudo-random pattern at high data rate. 



[Turn -on delay(S) 



PDF of turn on delay for pseudo random bit pattern at 1 


Gbi 


rn effect comes at high data rate. 


35 




■ J i It r s 1 u spc k cut ( bit r lk as o ig is t ic bit p lttc i is 010 
® Jitter increases due to bit pattern effect at high data rate. 

BS 1 he nature of variation of jitter with on-states is different for high data rate with pseudo 
random bit pattern, 

4 The practical relevance and the utility of the numerical 
mutation model: 


This simulation model approximately represents the real world system. 

It is important to investigate the turn-on jitter(due to spontaneous emission and bit pattern 
effect) since it causes the increase of BER at gigabit data rate in optical communication 
system. With this simulation model we can investigate the turn on behaviour of any single 
mode-single polarization VCSEL after introducing the parameters of the VCSEL in the 
program without performing the experiments. 

> 

The simulation results here matches well with the results obtained form the theoretical 
model described in [2], Since the simulation is performed by numerically solving the basic 
rate equations it actually verifies that analytical model is correct. 

From the simulation results we get how the jitter varies with the on-state level of the 
modulation current for a zero biased laser. So for high data rate to obtain less BER we can 
choose al what range of on-lcvel the laser should be operated to get BER less than 10' 9 
which is the tolerable BER for communication purpose. 

In simulation we can plot how the photon number and carrier number varies with the 
variation of current. Also we can observe the relaxation oscillation and we can investigate 
what parameters affects the oscillationfe.g gain compression coefficient k p (fig 4.4.1)). It 
can be observed that if k r is increased the oscillation decreases and for k p =25 W" 1 the 
oscillation is quite matching with the experimentally obtained oscillation. So we can say 
for this laser the value of k p -25 W' 1 approximately. 


In simulation by changing bit rate we can investigate the effect of bit rate on jitter and also 
we can measure upto what data rate the laser can operate because the upper limit of bit rate 
will depend on electron and photon turn-on delay also for zero biased laser. By applying 
high data rate( 1 Gbits/s) the effect of bit rate on jitter has been investigated. It shows that 
jitter does not depends on bit rate as long as 1010 bit pattern is applied. For pseudo randon 
the jitter changes with bit rate. Higher the bit rate high is the jitter. 

Here the simulation model has been done for source only. In optical communication systen 
we can simulate the transmission channel as well as the receiver aiso. Then one can 
connect different module by applying the output of one to the next and the whole 
communication system can simulated to investigate the behaviour of it. 

The eye diagram can be plotted by overlapping successive bits in this model and with the 
help of this diagram the BER can be measured for a particular data rate and for a particular 
on -stale. 





finalized Photon No.] 







® Also the effect of different bit pattern on jitter can be observed in this simulat on nn 
When the pseudo random pattern is app ied the simulation results shows that jitter 
increases significantly due to bit pattern effect especially for zero-biased operation 1 
effect has been verified experimentally also. Although the effect of bit pattern effect 
significant at very data rate( Gbits/s). 

B The simulation model can be modified and the work can be extended further for twt 
polarization multimode VCSEL. 
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Lab VIEW simulation programs 

1 . Main program (VCSELPRS2.vi) 

2. Parameter subprogram (parameter. vi) 

3. R-K subprogram (newRK4.vi) 

a) newfxy.vi 

b) newarr.vi 

c) newdSNdt dNdt.vi 

d) newGN.vi 

4. Noise subprogram (noi.vi) 

5. PDF subprogram (PDF-Direct. vi) 

6. Program for selection of parameters t e , t p (fr.vi) 
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